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ABSTRACT 

The evolution of the abundance of galaxy clusters depends sensitively on the value 
of the cosmological density parameter, f2 - Recent ASCA data are used to quantify 
this evolution as measured by the cluster X-ray temperature function. A \ 2 minimi- 
sation fit to the cumulative temperature function, as well as a maximum likelihood 
estimate (which requires additional assumptions about cluster luminosities), lead to 
the estimate flo ~ 0.45 ± 0.2 (1-er statistical error). Various systematic uncertainties 
are considered, none of which enhance significantly the probability that Qq = 1. These 
conclusions hold for models with or without a cosmological constant, i.e. with Ao = 
or 1 — Qq. The statistical uncertainties are at least as large as any of the individual 
systematic errors that have been considered here, suggesting that additional tempera- 
ture measurements of distant clusters will allow an improvement in this estimate. An 
alternative method that uses the highest redshift clusters to place an upper limit on 
Qq is also presented and tentatively applied, with the result that f2o — 1 can be ruled 
out at the 98 per cent confidence level. Whilst this method does not require a well- 
defined statistical sample of distant clusters, there are still modelling uncertainties 
that preclude a firmer conclusion at this time. 
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Galaxy clusters are the largest virialised objects in the 
universe and, as such, provide useful cosmological probes. 
For example, the evolution in the abundance of clusters is 
strongly dependent on, and is therefore a sensitive probe of, 
the cosmological density parameter, Qo, (e.g. Evrard 1989; 
Frenk et al. 1990; Lilje 1992; Oukbir & Blanchard 1992). 

The amount of observational data concerning high- 
redshift cluster properties has increased significantly in re- 
cent years. EMSS (Henry et al. 1992; Gioia & Luppino 1994), 
ASCA (e.g. Donahue 1996; Mushotzky & Scharf 1997; Henry 
1997) and ROSAT (Ebeling et al. 1997, Burke et al. 1997, 
Rosati et al. 1998) measurements of the X-ray emitting in- 
tracluster plasma now complement the low-redshift studies 
carried out by Edge et al. (1990) and David et al. (1993). 
In addition, the CNOC survey provides galaxy velocity dis- 
persions for a well-defined sample of high-redshift clusters 
(Carlberg et al. 1996) and further information on mass dis- 
tributions is coming from weak gravitational lensing mea- 
surements (e.g. Luppino & Kaiser 1997; Smail et al. 1997). 

The modelling of the cluster population has relied heav- 
ily upon the expression provided by Press & Schechter 
(1974) for the mass function of haloes. This has received 



support from dark matter simulations of structure formation 
(e.g. Efstathiou et al. 1988; White, Efstathiou & Frenk 1993, 
Lacey & Cole 1994). More recently the treatment of a bary- 
onic component in simulations has led to a closer link with 
observational measurements (Evrard 1990; Cen & Ostriker 
1994; Metzler & Evrard 1994; Bryan et al. 1994; Navarro, 
Frenk & White 1995, hereafter NFW). In addition, analyti- 
cal work investigating the evolution of properties of the clus- 
ter population as a whole has been performed by Kaiser 
(1986, 1991), Evrard & Henry (1991) and Bower (1997). 

There is now sufficient overlap between these strands 
of research to allow useful conclusions concerning cosmolog- 
ical parameters to be drawn from the available data (see 
however Colafrancesco, Mazzotta & Vittorio 1997). Several 
groups have already performed this operation (Henry 1997; 
Carlberg et al. 1997b; Bahcall, Fan & Cen 1997; Sadat, Blan- 
chard & Oukbir 1998; Blanchard & Bartlett 1997; Reichart 
et al. 1998) , although the lack of consensus on the measured 
value for f2o casts doubt on the reliability of at least some 
of these studies. The purpose of this paper is to investigate 
the systematic uncertainties in the estimation of fio using 
cluster evolution, and to understand why the various studies 
differ. 

In Section H the cosmology- dependent model for de- 
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scribing the evolution of the cluster population is described. 
The procedure for and results from fitting the models to the 
observed temperature functions are presented in Section ^| 
Section ^ contains a description of a maximum likelihood 
method used to estimate cosmological parameters using the 
cluster redshifts and temperatures. The results from these 
two methods are compared with other studies in Section |E| 
In Section ^ an alternative method is explored for using 
high-redshift cluster data to constrain Qo- A discussion of 
the different techniques is presented in Section [?]. The im- 
plications of these results for flo — 1 scenarios are discussed 
in Section ^, and conclusions are summarized in Section ^[ 



2 MODEL 

The model of the population of hot galaxy clusters employed 
here is the same one described by Eke, Cole & Frenk (1996). 
This is based upon the analytical expression for the comov- 
ing number density of dark matter haloes first derived by 
Press & Schechter (1974): 
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Here M v i r denotes the halo virial mass, defined such that 
a halo has a mean density A c times the critical density at 
redshift z, where A c is evaluated using the spherical col- 
lapse model; p represents the present-day mean density of 
the universe and a(M v i T ) the present-day, linear theory rms 
density fluctuation in spheres containing a mean mass M v „. 
Throughout this paper, we use the fit to the cold dark mat- 
ter power spectrum of Bardeen et al. (1986). The shape of 
the power spectrum and its amplitude at 8h _1 Mpc are de- 
scribed by as and V respectively (Efstathiou, Bond & White 
1992). The spherical collapse threshold, S c , is employed for 
most of the analysis in the rest of this paper. Bryan & Nor- 
man (1998) report that their work, combined with the ma- 
jority of published results, suggests that a slightly (< 10 per 
cent) higher value provides a better fit to mass functions 
from dark matter simulations. However, the recent simula- 
tions of Tozzi & Governato (1997) appear to support a lower 
value. This will be considered further in Section Monaco 
(1998) has provided a review of the present understanding 
of the mass function. 

An additional crucial ingredient in the model is the con- 
version from virial mass to emission-weighted X-ray temper- 
ature: 
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In this equation, /3tm is a conversion factor, the numerical 
value of which will be discussed shortly, and X is the hy- 
drogen mass fraction which is taken to be X — 0.76. If the 
specific galaxy kinetic energy equals the specific gas ther- 
mal energy within the virial radius, then a value of /3tm = 1 
would be appropriate for an isothermal gas. Alternatively, 
this value follows from applying hydrostatic equilibrium at 
the virial radius and assuming both that the sum of the 
logarithmic derivatives of the gas density and temperature 



profiles at this radius is —2 and that the gas temperature at 
this radius equals the overall emission-weighted value. Nei- 
ther of these reasons is particularly persuasive. Firstly, there 
is likely to be some residual kinetic energy in the g £LS elS 8l 
result of incomplete thermalisation (Evrard 1989). Secondly, 
numerical simulations show that the gas density is dropping 
off faster than r~ 2 at the virial radius and that the temper- 
ature is also decreasing rather than remaining constant (e.g. 
NFW). The value of /3tm chosen here is based solely upon 
the results of hydrodynamical cluster simulations. These 
usually include shock heating but no radiative cooling or 
heating from any source. Such experiments (NFW; Evrard, 
Metzler & Navarro 1996; Eke, Navarro & Frenk 1998; Bryan 
& Norman 1998; Frenk et al. 1998) suggest that /3 T m = 1.0 
is a reasonable choice, and it is adopted as the default value 
in this paper, although the dependence of the results on this 
assumption will be investigated in Section Q Values returned 
from different types of hydrodynamical simulations for a va- 
riety of cosmologies range between 1 and 1.3 (see Bryan & 
Norman 1998), with a slight preference for larger values of 
/3tm when Qo is larger. No significant mass dependence of 
/3tm is found. 

The redshift dependence of the mean value of /3tm for 
the 10 clusters simulated by Eke et al. (1998), and the scat- 
ter among the individual cluster values are shown in Fig. |l| 
Most importantly for estimating Qo from cluster evolution 
out to z ~ 0.4, there is no significant evolution in either 
the mean /3tm or in the magnitude of the scatter about this 
value. As the temperature function has such a steep depen- 
dence on cluster temperature, uncertainties and possible bi- 
ases in the mass-temperature conversion are very important 
(e.g. Pen 1998). A Gaussian scatter of 20 per cent in the 
temperature at a given mass has been included, consistent 
with Fig. |l], when calculating the model predictions. The ne- 
glect of physical processes associated with galaxy formation 
in hydrodynamical simulations of clusters is a concern, but 
the work of Metzler & Evrard (1994) suggests that even ex- 
treme feedback from stellar evolutionary effects cannot heat 
the intracluster medium of a rich cluster by more than about 
15 per cent. Similarly, NFW found that Ptm was decreased 
by about 20 per cent if preheating was included in rich clus- 
ter simulations in order to match the slope of the observed 
luminosity-temperature relation. The study by White (1991) 
comparing observed velocity dispersions with temperatures 
is also suggestive that non-gravitational processes have a rel- 
atively small effect on rich clusters. In Section [|, the possi- 
ble extent to which the results presented here might depend 
upon physics missing from the numerical simulations will be 
considered further. 



3 EVOLUTION OF THE CLUSTER 
TEMPERATURE FUNCTION 

3.1 Estimating the cumulative temperature 
functions 

The cluster X-ray emission-weighted temperature measure- 
ments that are used here come from two different sources. 
At low redshift, an updated version of the flux-limited 
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Figure 1. The top panel shows the evolution with redshift of the 
parameter Ptm> relating the emission weighted X-ray tempera- 
ture to the halo mass, as measured from an ensemble of 10 high- 
resolution hydrodynamical simulations of an Qo = 0.3, Ao = 0.7 
cosmology (Eke et al. 1998). The lower panel shows the standard 
deviation of the individual cluster Ptm values about this mean 
value as a function of redshift. 



(/(2 - 10 keV) > 3 x KT 14 Wm- 2 Q ) sample of 25 clus- 
ters compiled by Henry & Arnaud (1991, hereafter HA91) is 
employed (see Henry 1998 for details). These authors report 
that this sample is "at least 90 per cent complete". In this 
Section, a completeness of 100 per cent will be assumed, but 
the effect of this uncertainty will be investigated further in 
Section ^. The average ler measurement uncertainty in these 
temperatures is 7 per cent. For the high redshift data, the 10 
EMSS 0.3 < z < 0.4 cluster temperatures obtained by Henry 
(1997) are employed. These have an average la uncertainty 
of 14 per cent. No systematic differences are found between 
these results and those of Mushotzky & Scharf (1997) for 
the 5 cases where they studied the same clusters. 

In each case the cumulative cluster temperature func- 
tion is estimated using 



N(>kT) = J2 l l V ™ 



(3.1) 



kTi>kT 



where Vmax.s is the maximum volume in which cluster i could 

* For the benefit of older readers, 1 W = 10 7 erg s — 1 . 
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be detected given the survey selection criteria, and is calcu- 
lated using the cluster fluxes and redshifts. A complication 
for the high-z sample is that the area of sky surveyed was a 
function of the sensitivity, i.e. only a small fraction of the to- 
tal survey area was observed at the lowest flux levels. Using 
the data in table 3 of Henry et al. (1992), the following ex- 
pression has been found to give the area of sky, A, in which 
a cluster providing a detect cell flux, /dct, would have been 
detectable by the EMSS: 



1 - 3.05 e~ a41/dot + 2.30 e "°' 77/dct . 



(3.2) 



Here, fdet is in units of 10 _16 W m -2 and Atot is the maxi- 
mum area covered by the EMSS cluster survey, 735 square 
degrees. This fit is accurate to within 10 per cent at all de- 
tect cell fluxes above the Henry (1997) catalogue limit of 
2.5 x 10 _16 W m~ 2 . Using the cluster detect cell fluxes and 
redshifts provided in table 1 of Henry (1997), the cluster flux 
can be calculated as a function of redshift. U ma x,s is then 
evaluated by integrating the product of dV/dz and survey 
area at a particular redshift from z=0.3 to z=0.4. 

Before this estimate of N(> kT) can be compared 
with the model, a correction for incompleteness in the mea- 
sured high-redshift temperature function is required. This 
arises because of the non-zero scatter in the luminosity- 
temperature relation and the fact that the high-redshift sam- 
ple has a non-zero lower redshift limit. As a consequence, 
at the low redshift limit of the sample, the flux limit of 
2.5 x 10~ 16 W m -2 corresponds to a non-zero luminosity 
limit such that for any given temperature, only clusters with 
luminosities higher than this limit will be included. The 
incompleteness will be greater for cooler clusters because, 
on average, they have lower luminosities. To calculate the 
amplitude of this effect, the data of HA91 for low-redshift 
clusters were used to provide both a mean luminosity as 
a function of cluster temperature and a scatter about this 
mean. What observational evidence there is suggests that 
the luminosity-temperature relation evolves very little out 
to redshifts of 0.4 (Tsuru et al. 1996; Mushotzky & Scharf 
1997; Henry 1997). This gave 

0.3-3.5 
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3.541og 10 (fcT) - 2.53, 
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where Llj 7 is the 0.3 — 3.5 keV luminosity measured in 
10 37 h~ 2 W and kT is in units of keV, and a Gaussian scatter 
in logioL^ 3 " 3 - 5 of 



0.69 - 0.471og 10 (fcT). 



(3.4) 



At any given temperature, the fraction of clusters that are 
actually selected is then just given by 



F(T, Z/iim) = 0.5 erfc(a;ii m (T)) 



(3.5) 



where erfc represents the complementary error function, 
x'lim = (log 10 Ln m - log 10 Z)/(v / 2cr) and L llm is the 
(cosmology-dependent) luminosity of a cluster at z = 0.3 
having a flux equal to fu m . To allow for the fact that the 
catalogue flux limit refers to the EMSS detect cell whilst 
the clusters are typically more extended, a factor of 2.1 has 
been included to account for the conversion between these 
two fluxes (see figure 1 of Henry et al. 1992). The product of 
this factor and the differential number density at the same 
temperature can then be integrated with respect to temper- 
ature to give the correction to the cumulative temperature 
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Figure 2. The stepped curves show the low- (dashed) and high- 
(solid) redshift cumulative cluster temperature functions, calcu- 
lated assuming that Qq = 1. The points are the values that were 
used in the model fitting, and include the incompleteness correc- 
tion for the high-redshift case. Open and filled symbols correspond 
to the low- and high-redshift data respectively. Error bars were 
computed using bootstrap resampling. 

function. Even for temperatures as low as 4 keV the incom- 
pleteness is no more than about 20 per cent. This can be seen 
in Figure ^, where the low redshift cumulative cluster tem- 
perature function is shown as a stepped line along with that 
from the high-z sample, calculated assuming that flo — 1 
and ignoring the incompleteness correction. The points rep- 
resent the data used for the quantitative comparison with 
the models, and they include the correction factor which 
only has any significant effect on the lowest temperature 
measurement. The la uncertainties on the number densi- 
ties come from a bootstrap resampling procedure using 10 4 
catalogues. 

3.2 Comparison with model predictions 

Given that the cumulative cluster number density varies 
rapidly with temperature, uncertainties in either the mass- 
to-temperature conversion or in the observed temperatures 
of individual clusters will impact more on the results than 
any minor incompleteness. These effects can both be mod- 
elled by smoothing the model temperature function before 
it is compared with observational data. The width of the 
smoothing function, which has been taken to be a Gaus- 
sian, has been estimated by adding the fractional uncertain- 
ties from the two sources in quadrature. From the results in 
Fig. [jj, the scatter in individual cluster /3tm values is slightly 
less than 20 per cent and it appears to be independent of 
redshift for z < 0.4. While these simulations were of an 
Oo = 0.3, Ao = 0.7 model, in the absence of information 
to the contrary, this scatter is assumed to apply to other 
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Figure 3. The predicted evolution of the cluster temperature 
function between z = 0.05 and z = 0.33 is shown for an Oo = 1 
model with erg = 0.52 and T = 0.25, both before (solid) and 
after (dashed) a Gaussian smoothing with width 0.2T has been 
applied. 



cosmological models as well. For both observational redshift 
ranges, the temperature uncertainties are sufficiently small 
compared to the scatter in /3tm, that neglecting them will 
not affect the final results significantly. A Gaussian smooth- 
ing width of 0.2T was therefore used. From Fig. |^ it is appar- 
ent that this correction is not negligible compared with the 
amount of evolution that the Oo = 1 model predicts. How- 
ever, if the width of the smoothing Gaussian is similar for 
both the high- and low-redshift samples, then the amount of 
evolution is largely unaffected because the cumulative tem- 
perature function is almost a power law. The estimate of 
as will be decreased slightly by this effect. Compared with 
the rapid evolution in the number of clusters at a particular 
temperature, Fig. ^illustrates how little the temperature at 
a fixed number density evolves, and how important any evo- 
lution in the mass-temperature relation is over this range of 
redshifts. 

The best-fitting parameter set (Oo, I\ a$) was found 
by minimising \ 2 for the models calculated at five and three 
temperatures in the low- and high-redshift bins respectively. 
Bin temperatures were chosen to be far apart, in order to re- 
duce the correlations between the different points. Neverthe- 
less, the full covariance matrix was employed when calculat- 
ing the x 2 values. The number of bins used was the minimum 
for which the results of the model fitting were found to be 
robust to the exact positioning of the bins. If any fewer bins 
are taken, then the results depend upon their chosen tem- 
peratures. This is particularly so for the low-redshift data 
where the observed function is less smooth than in the high- 
redshift case. The statistical uncertainties in the best-fitting 
parameters were quantified using both the change in \ 2 an( A 
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Figure 4. The best-fitting model and the corresponding data 
are shown for the Ao = case. Low (high) redshift data are 
represented with a dashed (solid) line and open (filled) squares 
as in Fig. ^| 



the distribution of best-fitting values from the 10 4 bootstrap 
catalogues that were created from the observational data. 
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Figure 5. Histograms showing the distributions of the 10 4 boot- 
strap catalogue best-fitting Qo and as values. Open and flat mod- 
els are represented by solid and dotted lines respectively. 



3.3 Results 

The \ 2 minimisation gave rise to the following best-fitting 
parameter values: fio = 0.35, F = 0.07 and as = 0.71 
for A = and fi = 0.31, T = 0.07 and a$ = 0.78 for 
Ao = 1 — fio- In the latter case, two different effects act in 
opposite directions on the most likely value of fio compared 
to the Ao = case. If Ao is non-zero, then the expected 
amount of evolution at a particular fio increases because of 
the more rapidly changing growth factor. However, the vol- 
ume surveyed between redshifts 0.3 and 0.4 also increases 
if a Ao = 1 — fio term is included, and this has the effect 
of decreasing the observed cluster number density at high- 
redshift relative to the low-redshift measurement. Given that 
the most likely fio decreases, albeit only slightly, when Ao 
is assumed to be non-zero, the former effect is the more im- 
portant. The best-fitting curves for Ao = 0, together with 
the observational data points, are shown in Fig. W. 

In order to assess the statistical uncertainties in the es- 
timated parameters, one can consider both contours of A^ 2 
and the distributions of the best-fitting values of the boot- 
strap catalogues. It is interesting to look at both of these 
because the \ 2 method formally assumes that the uncer- 
tainty in the estimated temperature function is Gaussian 
distributed, which is not fully justified given the small num- 
ber of clusters in the two samples. Fig. ^ shows the distri- 
bution of best-fitting bootstrap catalogue parameters. The 
open and flat models are represented by solid and dotted 
lines respectively. Figures |andR show slices through the 



three-dimensional parameter space at different values of T 
and fio respectively. The contour levels represent 1, 2 and 
3<t contours for two degrees of freedom, i.e. not marginalised 
with respect to the parameter that is fixed. All panels are 
for open models with Ao = 0. Very little changes with the 
inclusion of a non-zero Ao. The minimum values of \ 2 are 
given for each of the slices. Estimates of the uncertainties 
have been obtained by looking at the cumulative distribu- 
tions of the bootstrap catalogue best-fitting parameters. The 
following ranges contain 68 per cent of these values: 

fio = [0.24, 0.49], T = [0.02, 0.16], as = [0.64, 0.76] (3.6) 

if Ao = 0, and 

fio = [0.20, 0.49], T = [0.02, 0.16], as = [0.67, 0.88] (3.7) 

if Ao = 1 — fio- When the individual bootstrap catalogue 
best-fitting parameters were overplotted on the A^ 2 con- 
tours shown in Figs. ^| and they were found to trace 
very similar patterns. Only a few lay outside the '3a' con- 
tours, having fio values slightly larger than would be ex- 
pected given these two figures. However, as is evident from 
Fig. |E[ there are nevertheless very few catalogues favour- 
ing fio > 0.9. This provides some reassurance that the A% 2 
contours are not providing a misleading impression of the 
statistical uncertainties in the various parameters, despite 
the low number of clusters per bin. 



© 0000 RAS, MNRAS 000, 000-000 



6 V.R. Eke et al 



1.5 - 



0.5 



1.5 



0.5 - 



1 1 1 1 

r=o.o2 : 

X z =3.8 - 

Amin 


1 1 1 1 1 1 II 1 II 1 1 

r=o.io : 

Y a . =2.4 - 

, Amin 

| i [ I | i | I I I | i j i | I | i 


T=0.25 I 
'. i v z - =5.3 

Amm 

[ 1 1 1 1 1 1 1 


Ii T=0.35 I 

— 1 — 

" \ Y a = 7 I 

- \ Amin 

1 1 ! 1 1 1 l 1 1 



0.5 



0.5 



Figure 6. Contours of constant Ax 2 for slices in parameter space 
at fixed r, as given in the individual panels. The A\ 2 values 
shown are 2.30, 6.17 and 11.8 and, for two interesting parameters, 
these correspond to enclosed probabilities of 1, 2 and 3<r. All 
panels are for the Ao = models. In each case, the minimum \ 2 
from all of the models with the chosen V is also given. 
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Figure 7. Contours of constant A\ 2 for slices in parameter space 
at fixed Co i as given in the individual panels. The A\ 2 values 
shown are 2.30, 6.17 and 11.8 and, for two interesting parameters, 
these correspond to enclosed probabilities of 1, 2 and 3<r. All 
panels are for the Ao = models. In each case, the minimum \ 2 
from all of the models with the chosen Qq is also given. 



4 MAXIMUM LIKELIHOOD PARAMETER 
ESTIMATION 

One drawback of the method presented in the previous Sec- 
tion is that the precise cluster redshifts are not employed and 
so the modelling ignores any evolution within the redshift 
ranges spanned by the two data sets. In this Section, details 
are given of a maximum likelihood technique that has been 
applied to determine the most likely model for producing 
the measured values of cluster temperatures and redshifts. 
This has the benefit over the previous method that no par- 
ticular binning in temperature is required and that the in- 
dividual cluster redshifts are used. The disadvantage is that 
the measured cluster fluxes are not included in the scheme, 
and additional assumptions about the cluster luminosities 
need to be made. 

If a(z,T)dzdT is the model probability of observing a 
cluster with temperature T at redshift z, then the likelihood 
of the observed distribution of cluster temperatures and red- 
shifts can be written as 



C = n^OidzdTe" 



n 



j=N+l c 



(4.1) 



Here j is an index running over all those of the N-rot to- 
tal bins in redshift and temperature that are empty, and i 
runs over just the N bins occupied by the observed clusters. 
This equation follows from assuming that the cluster red- 
shifts and temperatures are independent, and taking bins 
such that ai <C 1 Vi. Neglecting the terms that do no t de- 
pend upon the cosmological parameters, equation ( [l.l|) can 
be rewritten as 



ln£ 



L; =1 £li 



aAzAT. 



(4.2) 



(4.3) 



The function a(z,T) can be defined as 

a(z,T) = n{z,T,p) AV f ,P C(z,T,p) 
az 

where p represents the three model parameters being inves- 
tigated, n is the comoving cluster number density given by 
the Press-Schechter-based model including a correction for 
the 20 per cent uncertainty in the mass-temperature con- 
version, V is the comoving volume and ( is the product of 
the fraction of sky surveyed and the estimated completeness 
of the survey. For the high redshift catalogue, the fraction 
of sky surveyed has been calculated by integrating over all 
possible fluxes at a given (z,T,p), thereby accounting for 
the scatter in the luminosity-temperature relation and the 
variation of sky covered as a function of flux in the EMSS. 
A single temperature bremsstrahlung model has been used 
to calculate the 'K'-corrections to the cluster fluxes. The de- 
fault temperature to luminosity c onv ersions that have been 
employed are given by equation (3.3) for the 0.3 — 3.5 keV 



band and, again using the HA91 clusters, 



log io £3 



3.931og(fcT/keV) - 2.92 



(4.4) 



for the 2—10 keV band. The scatter in logio (luminosity) 
about these relations is described by equation (3.4) and by 



a = 0.70 - 0.521og 10 (fcT/keV) 



(4.5) 



respectively. To estimate how the choice of L — T conversion 
affects the results, the David et al. (1993) sample has also 
been considered, giving 
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Figure 8. A(— In likelihood)s for Co (left panel), T (middle) and <rg (right), marginalised over the other two parameters in each case. 
Dashed lines at 1 and 2a significance are shown, and the top of each panel corresponds to 3<r for one interesting parameter. The default 
assumptions giving rise to these results are described in the text. 



log 10 L^ 3 " 3 ' 5 = 2.761og(fcT/keV) - 2.02 (4.6) 
and 

logio^sf 10 = 3.451og(fcT/keV) - 2.57. (4.7) 

As in the previous Section, there is assumed to be no evolu- 
tion in the luminosity-temperature relation. 

The dependence of the best-fitting parameters on the 
following factors has been investigated: the assumed com- 
pleteness of the low-redshift sample, the assumed conver- 
sion between cluster flux and EMSS detect cell flux, the 
luminosity-temperature conversions, and the amplitudes of 
the scatter about these mean relations. The default choice is 
to assume that both high and low-redshift samples are 100 
per cent complete. Note that a lower completeness for the 
high-redshift sample would give rise to lower estimates for 
fio. Unless stated otherwise, a cluster at redshift z is taken to 
have 2.4 — z times as much flux as would fall into the EMSS 
detect cell (see figure 1 of Henry et al. 1992). The redshift 
and temperature limits for the binning are chosen so as just 
to encompass all of the data points. For the low-redshift sam- 
ple this means 2.5 <kT/keV< 10 and < z < 0.1, and for 
the high-redshift bins 3.7 <kT/keV< 10 and 0.3 < z < 0.4. 
The region of parameter space that has been searched for 
the most likely parameters is 0.01 < fio < 1 (with Ao = 
or 1 - fio), 0.01 < T < 0.25 and 0.5 < erg < 1.2. 

4.1 Results 

With the default assumptions described in the preceeding 
section, the most likely parameter values are fio = 0.52, 
r = 0.08 and erg = 0.63 for the Ao = model. The statistical 
uncertainties on these parameters can be determined by con- 
sidering the marginalised likelihoods. These one-dimensional 
likelihoods are shown for all three of the fitted parameters in 
Fig. 0. In each case the maximum likelihood is found at each 



value of the chosen parameter, independently of the other 
two parameters. Confidence levels are therefore taken for 
a single interesting parameter, leading to la uncertainties 
of fio = 0.52± ;^, T = OmtoZ and <r 8 = 0.63^;^. The 
model predictions for the marginalised number of clusters 
as a function of either redshift or temperature are shown in 
Fig. |^ for the most likely case. It is evident that the high- 
redshift sample cluster distribution is not very well matched 
by the model. In fact a Kolmogorov-Smirnov test indicates 
that the probability that this data set was produced from the 
model distribution is of order 10~ 3 . A comparison of CNOC 
(Carlberg, Yee & Ellingson 1997a) with EMSS (Henry et al. 
1992) redshifts suggests that the low probability cannot be 
accounted for by redshift measurement errors. If the maxi- 
mum redshift limit of the Henry (1997) cluster sample is de- 
creased from 0.4 to 0.33 and the likelihoods are recomputed 
then the resulting fio becomes <C 0.13 (the best-fitting erg 
is at least 1.2, which is the largest value considered on the 
grid of likelihoods). This suppression of the high-z survey 
volume with the corresponding removal of only one of the 
10 clusters has the effect of reducing the apparent evolution 
and, likewise fio. However, if the redshift distribution of the 
2 > 0.4 EMSS clusters is also considered, then models with 
so little evolution will overproduce distant objects so vastly 
that even the additional uncertainties in such predictions are 
unlikely to account for the results. The redshift distribution 
of the 0.3 < z < 0.4 clusters is therefore rather difficult to 
explain. 

Fig. [l(j shows the effect on the most likely parameters 
produced by altering the model assumptions. From top to 
bottom, the rows show the variation with 1) /3tm , 2) survey 
incompleteness and uncertain /ci us //det, 3) L — T treatment 
and 4) the inclusion of a non-zero Ao term. Table 2 contains 
details of the most likely parameters in the different cases. 

Increasing Ptm has the effect of boosting erg. This sup- 
presses the amount of evolution at any given fio , thus leading 
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Figure 9. Marginalised data (stepped curves) and best-fitting 
model (smooth curves) for the redshift and temperature distribu- 
tions of the clusters in the best-fitting model: Qo = 0.52 (Ao = 0), 
T = 0.08 and a s = 0.63. 



Table 1. Most likely parameters for the different input assump- 
tions shown in Fig. hoj 



Change from default choices 



H 



none 0.52 0.08 0.63 

f3 TM = 1.2 0.55 0.09 0.68 

f3 TM = 1.5 0.61 0.10 0.73 

(3tm = 2.0 0.72 0.11 0.81 

jW/det =2.0- z 0.59 0.07 0.61 

/cW/det = 2.8 - z 0.44 0.10 0.67 

2 = sample 90 per cent complete 0.58 0.08 0.62 

z = L - T scatter x0.7 0.59 0.11 0.61 

z = 0.3 - 0.4 L - T scatter x 1.3 0.56 0.07 0.62 

L-T given by eqns. (|J) and ( |0| ) 0.43 0.05 0.67 

z = 0.3 - 0.4 mean L\T X 1.5 0.64 0.06 0.60 

A = l-n 0.45 0.10 0.69 



to a higher estimate for the density parameter. It should be 
noted that a much more effective way of changing the most 
likely Qo would be to have some evolution in the value of 
/3tm, rather than merely changing the mean at all redshifts. 
Even /3tm = 2, which seems rather extreme (see Section^), 
only produces a +0.2 change in the density parameter. 



The second row in Fig. shows the relative insensitiv- 
ity of the results to the assumed conversion between total 
cluster flux and EMSS detect cell flux, Afi = ±0.08 for ex- 
treme choices (see Henry et al. 1992 figure 1). When the ratio 
/cius//dct is taken to be larger, this corresponds to a need 
for more luminous clusters in the high-redshift sample and 
hence a lower value for Qo- A 10 per cent incompleteness 
in the low-redshift cluster sample enhances the evolution 
enough to create Qo = 0.58. However, if the high-redshift 
sample was less than 100 per cent complete then the best- 
fitting value would decrease accordingly (see Section []). 

The scatter in the luminosity-temperature relations has 
been varied to see how much the most likely density pa- 
rameter can be increased. In the third row of Fig. hi it is 
apparent that very little impact is made with these alter- 
ations. More significant is the reduction in Qp th at co mes 
from exchanging the L — T e qua tions (4.6) and j4.7| ) for 
the default choices (|J) and (Q. This yields fio = 0.43. 
The amplification of the high-redshift luminosity at a given 
temperature by 50 per cent, which is more than is observed 
(Tsuru et al. 1996; Mushotzky & Scharf 1997; Henry 1997), 
increases the number of clusters that should be seen further 
away. Consequently, to fit the observed abundance, it is nec- 
essary to increase the amount of evolution to compensate, 
giving a best-fitting Oo = 0.64. 

As is evident from this figure, the best-fitting value of 
fio is difficult to change by more than ~ 0.1 with any sin- 
gle one of these systematic changes. Additional reasons why 
these results might be misleading will be considered in Sec- 
tion with particular emphasis on what can be said about 
the health of the critical density model in the light of cluster 
evolution. 

The final row in Fig. [l(j shows how the most likely pa- 
rameters are effected by the inclusion of a Ao = 1 — Qo term. 
In this case, the most likely parameters and la marginalised 
uncertainties are Q = 0.45±o re, r = 0.10l° o5 and °"s = 

69+ - 13 
u.oa_ og . 

Only the choice of upper redshift limit for the Henry 
(1997) sample and the excessive change in /3tm of the sys- 
tematic effects investigated above, manage to alter the most 
likely Qo by much more than 0.15. However, the statistical 
uncertainties are sufficiently large that in none of the above 
cases is Qo = 1 ruled out with a significance of 3cr or greater. 



5 COMPARISON WITH OTHER STUDIES 

The analysis presented here produces results in good agree- 
ment with those of Henry (1997). This is not surprising, 
given that very nearly identical data have been employed, 
albeit within different modelling procedures. These results 
fit happily with the ROSAT logN-logS analysis by Math- 
iesen & Evrard (1998) as long as the luminosity evolves like 
L oc (1 + z)~ 3 . Thi s is approximately what is implied by 
equations (2.2) and fe.3j), although it should be noted that 
are almost completely insensitive 
unlike those in 



the results in Section 

to the assumed form of the L-T relation, 
Section Q Kitayama & Suto (1997) have also used the logN- 
logS relation of X-ray clusters to arrive at the conclusion 
that Qo = 0.45 is favoured if Ao = 0. 

Of the other results in this field, there appear to be two 
more distinct camps, one on either side of the above settle- 
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Figure 10. A(— In likelihood)s for Qo (left panels), V (middle) and ag (right), marginalised over the other two parameters in each case. 
Dashed lines at 1 and 2cr significance are shown, and the top of each panel corresponds to 3<r for one interesting parameter. In all rows, 
the default curves shown in Fig. ^| are reproduced with solid lines. The top row shows curves for /3tm = 1-2, 1.5 and 2 with dotted, short 
dashed and long dashed lines respectively. In the second row, the results of a cluster to EMSS detect cell flux ratio of 2 — z (dotted) 
and 2.8 — z (short dashed) are given along with results obtained by assuming that the low-redshift data are 10 per cent incomplete 
(long dashed). Illustrated in the third row are curves resulting from decreasing the L — T scatter for the low-z sample by 30 per cent 
(dotted), i ncre asing t he s catter by 30 per cent for the high-z sample (short dashed), using the different mean L — T relations described in 
equations (4.7) and (4.6) (long dashed) or simply increasing the high-z luminosity at a given temperature by 50 per cent (dot-dashed). 
The final row shows the likelihoods found after assuming that Ao = 1 — f!o (dotted). 
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ment in the middle of f2o-space. On the lower-side, although 
not by much, favouring Q, « 0.40 ± 0.12 and 0.30 ± 0.10 
are Carlberg et al. (1997b) and Bahcall et al. (1997) respec- 
tively. These results are consistent with those found here, 
although the small quoted uncertainties strongly imply that 
fio < I- The 'high-ground', with Qo = 0.8 — 1, is occupied 
by Sadat et al. (1998) and Blanchard & Bartlett (1997). At 
this point, the reader might wonder why any of these con- 
clusions concerning flo are worth taking seriously, given that 
they appear to be inconsistent and span the whole range of 
interesting values. That such disparate results are obtained 
is presumably a suggestion that unaccounted-for systematic 
errors are lurking beneath the surface in at least one of these 
analyses. As has been stressed in this paper, and was also 
discussed by Pen (1998), the cluster number density evolu- 
tion is rapid if considered at a particular mass. However, at 
a fixed number density the mass varies relatively little over 
the range of redshifts probed by the presently available data. 
For instance, in an Qo = 1 model normalised to match the 
z — cluster abundance, the evolution between z = 0.05 and 
z = 0.33 of the mass at a constant comoving cluster num- 
ber density is ~ 45 per cent. Any unaccounted-for relative 
systematic temperature errors (i.e. affecting the low- and 
high-redshift samples differently) must therefore be smaller 
than ~ 30 per cent, or else the comparison between the 
models and observations will become completely unreliable. 
The quantification of any potential systematic uncertainties 
in mass estimation is thus crucial in determining the size of 
the error bar on the derived fio- 

For the studies producing lower-fio determinations, ve- 
locity dispersions obtained using different procedures are 
used to estimate masses. Frenk et al. (1996) and van Haar- 
lem, Frenk & White (1997) investigated the biases in mass 
estimation using velocity dispersions, and more recently 
Borgani et al. (1997) considered the model dependence of 
these biases. These studies suggest that only after very care- 
ful modelling of the velocity dispersion calculation proce- 
dure, should this method be applied with any confidence. 
From the simple scaling laws M oc cr 3 and M oc T 1 ' 5 , it 
is immediately apparent that fractional uncertainties in ve- 
locity dispersion would need to be half those in the X-ray 
temperature in order for the virial mass to be similarly un- 
certain. The importance of such conversions is shown by fig- 
ures 1 and 2 of Carlberg et al. (1997b) where the higher red- 
shift CNOC number density determination lies at a greater 
mass than the other CNOC point. This is slightly surprising 
because the virial theorem suggests that the same veloc- 
ity dispersion should be produced by equal masses out to 
a given physical radius, independent of redshift. Under the 
assumption that both CNOC redshift ranges have the same 
cluster mass limit, which was the choice made by Bahcall et 
al. (1997), the higher redshift CNOC point can be made to 
lie rather happily on the Qo = 1 model prediction shown in 
Carlberg et al.s figure 2, as long as their high-redshift point is 
shifted to the mass limit of their low-redshift sample. While 
it is not clear that this is the correct mass to choose, and 
Bahcall et al. took an intermediate value, this demonstrates 
how sensitively the conclusions about Qq depend on such 
details. In addition to these difficulties, the CNOC sample, 
which is almost a subsample of EMSS clusters and contains 5 
from the high-redshift set studied by Henry (1997), has both 
a luminosity and a velocity dispersion cut. Consequently, 



when converting to an observed number density of clusters 
above a particular mass, a correction for scatter in the cor- 
relations between these properties is required. As Carlberg 
et al. (1997b) showed, this is uncertain by ~ 50 per cent. 

Bahcall et al. also consider the richnesses in the Palo- 
mar Distant Cluster Survey. This is suspect because of the 
well-known difficulties associated with trying to calibrate the 
richness-mass correlation and, more importantly, any possi- 
ble evolution in it (e.g. van Haarlem, Frenk & White 1997). 

Moving on to the large-fio camp, Blanchard & Bartlett 
claim that their result "in conjunction with the absence of 
any negative evolution in the luminosity-temperature rela- 
tion, provides robust evidence in favour of a critical density 
universe" . The data used to draw this conclusion, come from 
a variety of sources: CNOC velocity dispersions, EMSS lu- 
minosities and ROSAT luminosities. As these authors point 
out, there is considerable uncertainty in the number den- 
sity of clusters above 6 keV at redshift zero. A consistent 
value of ~ 30 x 10 _8 (/i _1 Mpc)~ 3 is obtained from both 
the samples of Edge et al. (1990) and HA91. Integrating 
the ROSAT Brightest Cluster Sample (BCS) (Ebeling et al. 
1997) best-fitting Schechter function for the EMSS passband 
(0.3 — 3.5 keV) from infinity to the luminosity corresponding 
to 6 keV (1.25 x 10 37 /i~ 2 W) suggests a number density of 
clusters of ~ 40 x 10~ 8 (h^ 1 Mpc)~ 3 . The value calculated by 
Blanchard & Bartlett is twice as large as this. They derived 
their estimate by integrating in the ROSAT passband to a 
limit determined from the L-T relation of Arnaud & Evrard 
(in preparation). The consequence of having a larger red- 
shift zero abundance is to increase the apparent evolution 
and hence increase the derived fio- 

Sadat ct al. (1998) used the lack of evolution of the 
cluster luminosity-temperature relation to give a "tentative" 
measurement of Qo = 0.85 ± 0.2. This determination was 
based on an application of the method proposed by Oukbir 
& Blanchard (1997). They showed how, by using the redshift 
zero cluster temperature function normalisation and the red- 
shift distribution of EMSS clusters, it was possible to relate 
Q,o to the evolution of the cluster luminosity-temperature 
relationship. However, as Sadat et al. note, Oukbir & Blan- 
chard (1997) used the uncorrected (see Eke et al. 1996) 
HA91 temperature function, which gives an overabundance 
of present-day clusters. This will feed through into the f2o 
determination (equation 7 of Sadat et al.) and produce an 
overestimate, although it is not clear by how much. The 
measurement of the L-T evolution found by Sadat et al. re- 
lies upon the correct choice of the slope for this relation. 
It is important because the more distant observed clusters 
are intrinsically hotter and more luminous than the nearby 
ones, and if the slope of the L — T relation (they choose 
L oc T ) is not accurate then some spurious evolution (or 
lack of evolution) would be introduced because of the differ- 
ent types of clusters being probed at high and low redshifts. 
For the low-redshift clusters contained in their table 1 the 
mean temperature is 6.05 keV and the mean luminosity is 
15.9 x 10 37 W. From table 2 in their paper, the mean temper- 
ature and luminosity for clusters satisfying 0.4 < z < 0.6 are 
7.13 keV and 55x W 37 W. If the L-T relation is non-evolving 
then this suggests that L oc T~ 7 (This exponent is still > 6 
if the one high-redshift cooling flow cluster RXJ1347. 5-1145 
is removed from this calculation.), so there is an apparent 
inconsistency in the method that is likely to have an impact 
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on the derived value of fio- Their figure 3 and equation 5 
suggest that, if the high-z clusters are biased to have large 
luminosities, then this would cause fio to be overestimated. 



6 ORDER STATISTICS AND THE EMSS 
CLUSTER SAMPLE 

A number of recent papers have reported the existence of 
large mass concentrations at redshifts above 0.5 (Luppino 
& Gioia 1995, Donahue 1996, Luppino & Kaiser 1997, Don- 
ahue et al. 1998). Very often, these findings are accompanied 
by a claim that the existence of such objects proves trou- 
blesome for hierarchically clustering models of structure for- 
mation, particularly if Qo = 1- However, the magnitude of 
this difficulty is not quantified. The purpose of this Section 
is to propose and then apply a method by which the extent 
of the problem can be gauged. To this end, the question: 
'What is the expected redshift of the nth furthest galaxy 
cluster with a temperature above fcT?', will be addressed. 
The answer will depend sensitively upon Qo because of the 
dependence of the growth rate on flo- Thus, one should be 
able to use the observations of hot, distant clusters to dis- 
criminate between different values of flo- More specifically, 
given that any observational survey is likely to be less than 
100 per cent complete, the presence of a hot, massive cluster 
can be used to set an upper limit on fio- 



6.1 Method and application 

The Press-Schechter expression for the abundance of haloes 
as a function of mass has been tested using large N-body 
simulations out to about z = 0.5 for haloes having masses 
corresponding to fcT ~ 5 keV in an Qq = 1 universe (e.g. 
Eke et al. 1996) . It appears to provide a good description of 
the cluster mass function over this range of redshifts when 
the spherical collapse threshold density is used. The more 
recent, and larger, simulation of Tozzi & Governato (1997) 
suggests that the Press-Schechter equation overestimates the 
observed evolution in very massive haloes. However, it is not 
yet clear how much of this apparent discrepancy results from 
their particular choice of groupfinder. Whilst the theoretical 
expectation is still somewhat uncertain for very rare objects 
at z > 0.5, this is an issue that can be addressed in the 
near future with very large- volume dark matter simulations 
(e.g. Evrard et al. in preparation), so this need not be an 
insurmountable barrier to the usefulness of the approach 
described here. 

As was described in Section B the simple spherical col- 
lapse model will be used in conjunction with the Press- 
Schechter expression in order to provide a description of 
the halo mass function. Additionally, /3tm = 1 will be as- 
sumed to hold at all redshifts, along with a 20 per cent scat- 
ter about this mean value. Using the results of Eke et al. 
(1996), the normalisation of the mass fluctuation spectrum 
will be taken as 0.50 ± 0.04 for the fio = 1 model considered 
in this Section. (Note that this is slightly different from the 
best-fitting value of 0.52 found by Eke et al. (1996). This 
new estimate takes into account both the systematic over- 
estimation of as resulting from the measured temperature 
uncertainties (see table 1 of Eke et al.) and another, slightly 



larger, term (neglected in their table) arising from scatter in 
the mass-temperature relation.) 

From these stipulations, it is possible to calculate the 
redshift distribution of all clusters hotter than a specified 
value. The probability that, from the population of clus- 
ters hotter than a particular temperature, a cluster is sit- 
uated at a redshift no greater than z will be denoted by 
P(< z) = P(< z\kT > fcT min ). This quantity is calculable 
from the model for any given temperature. For brevity and 
clarity the temperature dependence is suppressed in the fol- 
lowing expressions. Using this notation the probability that 
the furthest cluster (at z\) has a redshift greater than z can 
be written 



P(zi > z) = 1-P(zi < z) 



(P(< z)Y 



(6.1) 



where N is the total number of clusters above the cho- 
sen temperature limit in the whole volume, given by the 
Press-Schechter formula. For a particular redshift distribu- 
tion, larger N will probe further into the tail of the distribu- 
tion and zi will be expected to be larger. This procedure can 
be extended to give the redshift distributions of the fcth fur- 
thest cluster in the following manner. The probability that 
the fcth furthest cluster (hotter than fcT m i n ) has a redshift 
greater than a particular value is given by 



P(z k >z) = l- P{z k < z) 



(6.2) 



To calculate the last term in this equation, it is necessary 
to consider the various combinations of clusters that would 
give rise to z k < z. This probability can then be written, 
using the combination function C Q as 



P(z k < z) 



\ ~* jv. 



d P(< zf- j (l-P(< z)Y 



(6.3) 



and this, together with equation (6.2) gives P(zk > fc). This 
probability will be very sensitive to both the model normal- 
isation (78 and the error in the measured cluster temper- 
atures. Consequently, a marginalisation technique allowing 
for Gaussian errors on both of these quantities has been em- 
ployed. This involves integrating over all normalisations and 
temperatures, weighting with the probabilities that these are 
in fact the actual values of as and the lowest temperature 
of the clusters under consideration. 

The three clusters that will be used in this analysis are 
M50015.9+1609 at z = 0.54, M5"0451.6- 0305 at z = 0.54 
and A'/S1054.5 - 0321 at z = 0.83. ASCA temperatures of 
8.0 ± 0.6 keV (Mushotzky & Scharf 1997), 10.4 ± 1.2 keV 
(Donahue 1996) and 14.7 ± 2.1 keV (Donahue et al. 1998) 
respectively have been found for these objects. Therefore in 
the EMSS cluster sample, the third furthest cluster having 
a temperature greater than or equal to 8.0 ± 0.6 keV has a 
redshift larger than 0.53. Given that these clusters were all 
found during a search containing only 735 square degrees of 
sky, this value, rather than the entire 4-7T steradians will be 
used for calculating the model predictions. 



t The combination function is N Oj = N\/(j!(N - 
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Table 2. The main strengths and weaknesses of the three methods used here for constraining Qq using cluster evolution. 



Method Strengths Weaknesses 

X 2 (Section ^) Utilises the cluster flux measurements. Requires carefully selected data. 

Almost independent of L — T relation. Results depend slightly on z and T bin choice. 

M.L. (Section ^) All cluster redshifts are used. Requires carefully selected data. 

No z or T bin choices are necessary. Uncertainty introduced with L — T relation. 

Order statistics (Section ^) No need for a statistically well-defined data sample. Based on a model extrapolation. 

Independent of L — T relation. Results are very sensitive to the assumed erg. 
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Figure 11. The probability that the fcth furthest of all EMSS 
clusters with a temperature exceeding 8.0±0.6 keV has a redshift 
larger than z, assuming a cosmological model with f2rj = 1, erg = 
0.50 ± 0.04 and V = 0.25. There is less than a 1 in 7 chance that 
three such objects will exist in an EMSS cluster survey area of 
sky, and only a probability of ~ 0.02 that the third furthest hot 
cluster will have a redshift greater than 0.53. 



6.2 Results 

As can be seen in Figure [H], the probability that three 
clusters with kT > 8.0 ± 0.6 keV exist in an Qo = 1, 
erg = 0.50 ± 0.04, T = 0.25 model is less than 0.14. Fur- 
thermore, it is another factor of ~ 6 less likely that the 
third most distant of these has a redshift of more than 0.53, 
as is observed. The conclusion from this exercise would ap- 
pear to be that such a model can be ruled out at the 98 per 
cent confidence level. This would improve to ~ 99 per cent 
if an extra EMSS cluster hotter than 8 keV were found at a 
redshift larger than 0.53. It is worth reiterating that this re- 
sult would vary somewhat if the Press-Schechter treatment 
of the halo mass function was badly awry, or if the best- 
fitting erg were underestimated because of a poor choice of 
the factor relating mass to temperature, /3tm- Therefore, 
it would be premature to consider the existence of these 
hot clusters at high redshifts as being sufficient to rule out 
Qo = 1 cosmologies. Performing the above analysis for an 



Qo = 0.52 and erg = 0.63 ± 0.06 model (keeping F fixed at 
0.25) yielded a probability of 0.38 for the third furthest of 
the EMSS kT > 8.0 ± 0.6 keV clusters having z > 0.53. 



7 OVERVIEW OF THE DIFFERENT 
TECHNIQUES 

Having analysed cluster data in three distinct ways, it is 
worth discussing the strengths and weaknesses of the dif- 
ferent methods that were employed in Sections ^, ^ and ^. 
While all of these techniques use similar input data in an at- 
tempt to determine Qo, they are sufficiently diverse to merit 
some discussion. Table 1 contains a summary of the main 
pros and cons of these approaches. 

The main drawback with the first two methods, which 
is not so important for the order statistics approach, is the 
need for a statistically well-defined sample of data with well- 
understood selection effects. Given such observational cat- 
alogues, the particular data that are used by the \ 2 an( A 
maximum likelihood (M.L.) techniques differ somewhat. In 
the former case, only the cluster temperatures and fluxes 
are important. The 1/V ma x estimator enables this method to 
sidestep the need for any cluster luminosity-temperature re- 
lation, apart from a minor incompleteness correction. How- 
ever, the redundancy of the individual redshift measure- 
ments, other than to determine if a cluster should be in- 
cluded in the 'high-' or 'low-'z temperature function, leaves 
a certain freedom to choose what redshifts to assign to the 
fits. The choice of temperature bins is also something that 
can alter the resulting Qo determination to a small extent. 
In contrast, the M.L. method uses all the cluster tempera- 
tures and redshifts and avoids awkward binning problems. 
However, the limits of integration in equation (4.2), coming 



from the adopted ranges in z and T, do create some freedom 
in the most likely parameters. This comes as a result both of 
the peculiar redshift distribution of the high-z clusters, and 
the uncertainties in the selection functions at the extreme 
temperatures and redshifts. Additionally, the price paid for 
not using the measured flux values is that one must adopt 
a luminosity-temperature relation (at all relevant redshifts) 
and any error in the assumed form of this relation will in- 
troduce a systematic error into the estimated Qo. 

The difference between the derived best-fitting param- 
eters from the \ 2 an< A M.L. techniques has non-negligible 
contributions from the systematic uncertainties mentioned 
in the previous paragraph. Looking at the situation more 
positively, it is perhaps reassuring that, despite the various 



© 0000 RAS, MNRAS 000, 000-000 



Q from cluster evolution 13 



differences between the two methods, such similar results 
were returned. 

An opportunity to win the jackpot is offered by the or- 
der statistics method; namely the chance to rule out Qq — 1 
with observations of only a few very hot and distant clusters. 
However, the size of the payout for an apparently successful 
data set is not yet well-defined because of the uncertainties 
in both the cluster mass function evolution and the mass- 
temperature conversion for such rare objects at such high 
redshifts. Nevertheless, these are both eminently measure- 
able using numerical simulations, and the Hubble volume 
simulation performed by the Virgo consortium (Evrard et 
al. in preparation) should address the former of these issues 
in the near future. It is worth noting that the potential non- 
virialisation of large mass concentrations at high redshifts, or 
the boosting in emission- weighted temperature coming from 
a recent merger event should not cause any difficulties for 
this method. This is because such effects are included in hy- 
drodynamical simulations and would be calibrated through 
the mean /3tm and the scatter around this value. 



8 WHAT WENT WRONG? 

Given that all 'right-minded' cosmologists know that flo = 1, 
a pertinent question to ask at this point, on their behalf, is 
"What has gone wrong?", i.e. what would it need for an 
flo = 1 model to be more favoured by cluster evolution. 
In this Section, consideration is given to the magnitudes 
of the errors that are required in some of the preceeding 
assumptions. 

If flo is in fact unity, then it is very probable that at 
least one of the observations or the model used above is 
wrong. The first of these possibilities seems unlikely. An 
flo = 1 model would require a systematic overestimation 
of the high-redshift cluster temperatures by ASCA of about 
25 per cent, or a corresponding underestimation of the low- 
redshift temperatures by this same amount. This corre- 
sponds to a shift in all of the cluster temperatures that is 
typically twice the quoted uncertainty and, as such, seems 
highly unlikely. Having said that, the analysis of ASCA data 
by Markevitch et al. (1998) showed, for single temperature 
fits to 15 of the original 25 HA91 clusters, an average 6 
per cent increase in the measured temperatures. While this 
is not sufficiently large, if such a systematic shift were ap- 
plied to the entire low-redshift temperature function, then 
the most likely flo would increase by ~ 20 per cent. 

Contamination of the X-ray emission by active galac- 
tic nuclei (AGN) or cooling flows would alter the mass- 
temperature conversion. The effect of AGN contamination 
is to increase the emission-weighted temperature. Cooling 
flows have the opposite effect. These potential complications 
are better addressed from the point of view of decontaminat- 
ing the data, rather than adding extra processes into the nu- 
merical simulations, so will be considered as 'observational' 
problems. The likelihood of any interloping AGN in the high- 
redshift sample is very low, and high-resolution images show 
no evidence for such contamination. More important is the 
effect of cooling flows on the two different redshift samples. If 
cooling flows were more prevalent in the low-redshift sample 
than the high-redshift one, then their removal would result 
in an increased separation of the two temperature functions, 



thus favouring higher flo - By fitting both single-temperature 
and two isothermal models to ASCA data, Markevitch et al. 
(1998) showed that, for the 15 low-z clusters in common with 
HA91, the decontaminated emission-weighted temperature 
is, on average, 6 per cent larger than the single-temperature 
value. For the most discrepant single cluster, the effect was 
5 times as great. Allen (1998) performs a similar analysis 
including two of the Henry (1997) (i.e. high-z) clusters, and 
finds the temperatures increase by 15 and 24 per cent when 
a cool component is removed. There are two simple reasons 
why one might expect cooling flows to be more prominent 
in the EMSS high-z sample rather than the HA91 catalogue. 
Firstly the densities are larger and hence cooling times are 
shorter in the higher redshift range. Secondly, the EMSS se- 
lection algorithm tends to favour clusters that have centrally 
peaked emission, with only those exceeding a flux threshold 
in a 2'. 4 x 2'. 4 detect cell being included in the catalogue. If 
the results from Allen (1998) are typical of the whole Henry 
(1997) sample, then the effect of cooling flow decontamina- 
tion would be to increase the high-redshift cluster tempera- 
tures by more than the low-redshift ones. This would some- 
what overcompensate for the systematic temperature shift 
described in the previous paragraph, leading to an overall 
estimate of flo that is about 20 per cent lower than the 
value found in Section ^, where these considerations were 
neglected. 

As was discussed in Section ^ a 10 per cent incom- 
pleteness in the low-redshift sample only increases flo by 
0.06, so this is not a sufficiently large effect to help the crit- 
ical density model. Independent checks are available on the 
sample completenesses. These come from the comparison of 
luminosity functions obtained from the cluster catalogues 
used here with those from the larger ROSAT cluster sam- 
ples at low (Ebeling et al. 1997) and high (Rosati et al. 
1998) redshifts. ROSAT does not detect enough very lumi- 
nous clusters to compare with the majority of the luminos- 
ity range of the EMSS high-z sample, because of the smaller 
fraction of sky it surveyed. Nevertheless, where an overlap 
in luminosity exists, there is agreement between the clus- 
ter abundances in the different surveys. These comparisons 
have enough statistical uncertainty to make them reassur- 
ing rather than highly restrictive. There could still plausibly 
be ~ 40 per cent difference between the various samples. In 
order to favour flo = 1 though, there would need to be a sig- 
nificantly greater incompleteness in the low-redshift sample 
than exists in the more distant one, and this seems unlikely. 

The considerations in the preceeding few paragraphs 
suggest that if anything, relative to the results found in Sec- 
tions U and ^j, lower rather than higher values of flo would 
be favoured by the potential 'observational' systematic ef- 
fects. Thus, if the data used in this paper are the product 
of an flo = 1 universe then it seems that the finger of suspi- 
cion should point firmly in the direction of the modelling of 
cluster evolution. 

While the statement that the Press-Schechter expres- 
sion provides a good description of the mass function of 
galaxy-cluster-sized objects had stood up to much scrutiny, 
the recent report by Tozzi & Governato (1997) suggests 
that it overestimates the actual evolution for flo = 1- Their 
simulation is of a cube with side 500/i _1 Mpc and thus al- 
lows the evolution of very rare objects to be traced back 
to higher redshifts than was possible with the smaller vol- 
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Figure 12. A(— In likelihood) for Uq (left panel), T (middle) and 
erg (right), marginalised over the other two parameters in each 
case. Dashed lines at 1 and 2<r significance are shown, and the 
top of each panel corresponds to 3<r for one interesting parameter. 
The curves are calculated assuming that the collapse threshold is 
given by 5 c ff = 1.48(1 + z)~ ' 06 rather than the spherical collapse 
<5c 



umes used in previous studies. These authors suggest that 
the spherical collapse threshold S c should be replaced by 
S a ff = 1.48(l + z)~ 0,06 . This result appears to be slightly dis- 
crepant with previous findings where comparison is possible, 
but it may be that this is just an artefact of the groupfinder 
that was employed (see Cole & Lacey 1996 for more details) . 
Nevertheless, it is worthwhile to investigate the possibility 
that an evolving threshold is required. Assuming that 8 e ff 
is the form of S c that is applicable to all Ao = models 
(similar redshift evolution is found in a large Qo = 0.3 CDM 
simulation; Governato, private communication), the effect of 
such a varying threshold on the most likely value of Qq can 
be investigated using the likelihood procedure described in 
Section EL The resulting marginalised likelihoods are shown 
in Fig. [12l It is apparent that, whilst the most likely value of 
Qo has only increased to 0.58, the A(— ln£) required to reach 
fio = 1 is now just 2.2, i.e. only ~ 2a. Whilst this gives the 
impression that a small hit has been scored on the good ship 
fio < 1, the captain might retaliate by pointing out that the 
EMSS survey found clusters with z > 0.4 and, given that 
the Tozzi & Governato 8 e s fitted their simulated mass func- 
tion out to a redshift of one, these distant objects can now 
be included in the likelihood calculation. Temperatures have 
yet to be measured for all of these clusters, but they do have 
known redshifts, so an additional contribution to the like- 
lihood is computed by replacing the a in equation (4.2) by 
the integral of a over all cluster temperatures (2 — 20 keV 
in practice) and just comparing the probability of finding 



clusters at the observed redshifts. Additionally, the value of 
/cius//dct is taken to be 1.75, independent of redshift for the 
z > 0.4 EMSS clusters. The 90 identified EMSS clusters with 
2 > 0.4 and their corresponding redshifts are taken from 
Gioia & Luppino (1994) and Luppino & Gioia (1995), and 
the marginalised likelihoods resulting from their inclusion 
are shown in Fig. [H| It should be noted that there are also 9 
unidentified EMSS sources (Gioia & Luppino 1994), some of 
which might contribute to this list, and consequently reduce 
the estimated value of £Iq from that found here. The most 
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Figure 13. As for Fig. |12| but with the redshift information for 
the 9 z > 0.4 EMSS clusters included in the likelihood calculation. 



likely parameters change only slightly, becoming Qo — 0.62, 
r = 0.06 and as = 0.53. Furthermore, the inclusion of the 
more distant clusters significantly decreases the sizes of the 
statistical uncertainties on the most likely parameters, such 
that fio = 1 is now excluded at a higher confidence level than 
was the case with the original default assumptions, despite 
the favourably evolving 8 c g. Distributions of the numbers of 
clusters marginalised with respect to either redshift or tem- 
perature are shown in Fig. [w]. As the temperatures have yet 
to be measured for the highest redshift clusters, this model 
'prediction' is shown as a dotted line. 

In addition to any uncertainties in the accuracy of the 
Press-Schechter description of the mass function, there is 
some leeway in the conversion from mass to temperature 
that may effect the most likely Qq- There are at least 5 
options for tinkering in an attempt to salvage £lo = 1: 

1. Evolution of the scatter in both the cluster 
mass-temperature conversion and the observational 
temperatures would give rise to a systematic error in the 
estimation of Qo. It has been assumed that the intrinsic 
scatter is unchanging and the low- and high-redshift 
temperature errors are similar (which the quoted values 
suggest they are) . If these scatters at high redshift implied 
that a Gaussian smoothing width of 0.4T rather than 0.2T 
was appropriate, then even the Slo = 1 model cluster 
temperature function would appear not to evolve. 
However, for the mass-temperature relation, the low-S7o 
simulations of Eke et al. (1998) show no significant 
evolution in the scatter over this range of redshifts (see 
Fig. |l]). The observational uncertainties are relatively small 
and, only slightly larger for the high-z sample. The 
conclusion is therefore that such an escape for Qo = 1 is 
not easily arranged. 

2. Evolution of fj,, the mean molecular weight of the 
intracluster gas, has been ignored in the modelling. If it 
were 25 per cent larger at z = 0.33 than z = then the 
model predictions for fio = 1 would show no evolution over 
this redshift range. This change in \i is both extreme and 
in a counter-intuitive direction. Mushotzky & Loewenstein 
(1997) recently used ASCA data to show that the iron 
abundance in the intracluster medium (ICM) does not 
evolve over the range < z < 0.3. Also, Kauffmann & 
Chariot (1998) adopted a semi-analytical approach to the 
formation and evolution of elliptical galaxies, concluding 
that the metallicity of the ICM is approximately constant 
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Figure 15. The temperature functions for an Qq = 1, ag = 0.52 
and r = 0.25 model are shown with solid and dashed lines for 
z = 0.05 and z = 0.33 respectively. In each case, the curve s re - 



sulting from applying the temperature mapping of equation (S.l) 



are also shown. These are steeper than the original predictions. 
Stepped curves show the observationally determined temperature 
functions. 



Figure 14. Marginalised distributions of clusters. Calculated as- 
suming that the threshold varies with redshift according to 8 e g , 
and using the 9 z > 0.4 extra EMSS clusters. Smooth curves are 
shown for the most likely model, namely Qo = 0.62, T = 0.06 and 
erg = 0.53. The dotted curve in the temperature panel corresponds 
to the model prediction for the distribution of the temperatures 
that are yet to be measured. The stepped curves represent the 
observed cluster distributions. 



out to z > 1. This approach for saving Qo — 1 can be ruled 
out. 

3. Evolution of /3tm to the extent described for /i in point 
2, would also make Qo — 1 compatible with the 
observations. The definition of /3tm is such that it would 
need to be 25 per cent smaller at the larger redshift. As 
shown by the simulations of Eke et al. (1998) (see Fig. |l|) 
this amount of evolution is difficult to justify and, if 
anything, the (3tm at z — 0.33 is a few per cent larger than 
that at low redshift. For the various cosmologies studied 
numerically by Bryan & Norman (1998), a similar lack of 
evolution was also found (see their figure 4). This suggests 
that such a route will not assist Qo = 1 models either, and 
that the most likely value for Qo quoted here may be 
slightly biased to a value that is too large. 

4. One clear deficiency of the hydrodynamical simulations 
that have been used to relate virial mass to 

emission- weighted X-ray temperature is the lack of any 
heating of the intracluster gas by supernovae. That this is 
a problem becomes apparent when the simulated cluster 
luminosity-temperature relation is compared with the 



observed one. Whilst the models and scaling laws show 
L oc T 2 , the observations suggest a steepe r de pendence of 
luminosity on temperature (see equation (3.3)). Allen & 



Fabian (1998) find that decontaminating the cooling flow 
contribution to the X-ray emission can leave L oc T~ 2 . 
However, a similar study by Markevitch (1998) implies 
that, in order to produce this shallow slope, some feedback 
is still necessary after the cooling flow correction has been 
applied. NFW showed that by preheating gas in the 
manner suggested by Evrard & Henry (1991), the 
low-temperature simulated clusters could be warmed up 
enough, along with the corresponding decrease in 
luminosity as a result of the less dense cores, to produce a 
sufficiently steep luminosity-temperature relation to match 
observations. To obtain a very crude idea of the maximum 
extent to which such a change will affect the evolution of 
the Qo — 1 temperature functions, a simple mapping of 
temperature has been applied to the no-preheating value, 
(fcT)o, to give a 'feedback' temperature, (fcT)i. The 
mapping that has been adopted is 

log 10 (fcT)i = 0.571og 10 (fcT) + 0.34. (8.1) 

This route for converting the simulated clusters to lie on 
the observed luminosity-temperature relation is extreme in 
its effect on the cluster temperatures. In practice, one 
would expect any energy injection to heat the gas and 
decrease the central density. As a result, the simulated 
clusters would fall onto the observed L x — T x relation more 
as a result of decreasing luminosity than increasing 
temperature as is assumed for the mapping above. 
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Figure |l^ shows the predictions before and after 'feedback' 
when the same distortion is applied to both high- and 
low-redshift models. With the 'feedback' essentially just 
heating the cooler clusters, it produces steeper cumulative 
temperature functions, thus favouring lower V values. (It is 
worth noting that, if the lower temperature clusters were 
preferentially contaminated by cooling flows then the 
estimated F would be biased low, but Markevitch (1998) 
reports that slightly the reverse is seen in his sample of 
nearby clusters.) Interestingly, the extreme feedback 
simulations of Metzler & Evrard (1994) increased the gas 
temperature of a 5.6 keV cluster by ~ 15 per cent. This 
size of fractional change is only reached with the above 



mappin g for clusters cooler than ~ 4.5 keV. Despite all of 



this, the predicted amount of evolution still exceeds that 
observed. Only if such a prescription were proscribed for 
the high-redshift sample alone would feedback be effective 



in redu cing the amount of evolution seen in flo — 1 models. 
5. The bffocto of having a two phaoo intracluator gaa 



medium or magnetic fields are also lacking from the 
hydrodynamical simulations used to normalise the 
mass-temperature relation. While it seems unlikely that 
they should have a different effect on the two redshift 
ranges considered here, they could still bias the best-fitting 



parame ters as was illustrated in Section i , by shifting a 



different part of the mass function into the observed range 
of temperatures. The amount by which /3tm can be 
affected is roughly constrained by comparisons between 
weak lensing masses and X-ray masses. Present results 
(Smail et al. 1997; Allen 1998) appear to imply that this 
additional physics is not sufficiently misleading to allow 
flo = 1 to be consistent with the lack of evolution of the 
cluster number density. 



9 CONCLUSIONS 

In this paper, the value of the density parameter flo has been 
constrained using three different, yet similar, methods. The 
results are entirely consistent with those presented by Henry 
(1997). Combining the results from Sections ^ and |i| gives 
the most likely parameters as fl « 0.44±0.2, T « 0.08±0.07 
and 



as |~ 0.07 ±0.1 if Ao — 0. The 1 inclusion of a noii-zeio 
cosmolo gical conatant term ha3 very little effect, leading to 
flo « 0.38 ± 0.2, T « 0.09 ± 0.08 and cr 8 » 0.74 ± 0.1. 
Apart from the difficulty arising because of the redshift dis- 
tribution of the Henry (1997) sample, none of the systematic 
effects that have been considered alter the best-fitting flo by 
more than one of the statistical standard deviations quoted 
above. Unless an important systematic error or uncertainty 
has been overlooked in the analysis, these results would ap- 
pear to be limited by a lack of data, so, while this prelimi- 
nary study suggests that the flo = 1 model is still bobbing 
about unhappily in the water at a rejection level of 2 — 3a, 
a larger sample of clusters would offer the opportunity, with 
this method alone, to sink it properly. Enough ammunition 
should be provided by the upcoming X-ray telescope mis- 
sions XMM and AXAF. 
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